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摘要 :为 研究 锚 杆 拉 拔 实验 中 锚 杆 与 岩 体 间 的 力学 作用 过 程 和 机 制 ,在 非 连续 变形 分 析 (discon- 
tinuous deformation analysis ,DDA ) 方 法 中 实现 锚 杆 - 岩 体 界面 耦合 作用 算法 ,采用 多 线性 和 妹 结 滑 移 界 
面 模型 模拟 锚 杆 的 滑 移 脱 黏 ,并 针对 锚 杆 拉 拔 力作 用 节点 和 脱 黏 节点 提出 专门 算法 。 Es 
X5 SC 2E Rd xp vo sp sc dk 85 m ACER va ALPE SE E TRUST d TO E 4 CIE IE L9 HUE PT S 
HEUS GER Hkh E EE R E, i8 b CTAAEGEAA EA $ R EGE 2k i SNC S 


IPAR DDA 模拟 ,研究 了 锚 杆 的 拉 拔 力 - 拉 拔 位 移 关系 ,以 及 不 同 拉 拔 力 条 件 下 锚 杆 - 岩 体 界面 切 


) 应 力 和 锚 杆 轴 力 的 分 布 与 变化 规律 ,并 探讨 了 锚 杆 弹性 模 量 对 锚 杆 拉 拔 力学 响应 的 影响 。 锚 杆 - 岩 
s 体 界面 者 合作 用 算法 的 实现 ,也 为 岩 体 锚固 的 DDA 模拟 分 析 创造 了 条 件 。 

关键 词 : 锚 杆 拉 拔 ; 黏 结 滑 移 模型 ; 轴 力 ;界面 切 应 力 ; 非 连 续 变形 分 析 
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Abstract : To study the mechanical interaction process and mechanism between rock bolt and rock mass in 
the pullout experiment , an interface coupling algorithm between the rock bolt and rock mass is introduced 
into the discontinuous deformation analysis ( DDA) method , and the slip debonding of the rock bolt in the 
rock mass is represented by multi-linear bond-slip interface models. By comparing the simulation results 
with the experimental results, the parameters of the tri-linear and four-linear bond slip interface models 
from literatures are revised , and the correctness and applicability of the proposed interface coupling algo- 
rithm are validated. Through DDA simulations of rock bolt pullout based on the multi-linear bond-slip in- 
terface models with revised parameters, the pulling force and displacement relationship , and the distribu- 


tion and variation rules of the interface shear stress and bolt axial force under different pulling force levels 
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are studied. The influence of the rock bolt modulus on the rock bolt pullout mechanical response is also in- 


vestigated. 


Key words :rock bolt pullout; bond-slip model; axial force; interface shear stress; discontinuous deforma- 


tion analysis 


岩 体 的 错 杆 加 固 广 泛 应 用 于 地 下 开 挖 工程 、 岩 
体 边 坡 工程 和 矿山 工程 等 2 。 岩 体 锚 固 通过 锚 杆 
将 载荷 从 不 稳定 岩 体 区 传递 至 稳定 岩 体 区 ,力学 作 
用 过 程 和 机 制 复 杂 , 涉 及 锚 杆 与 岩 体 ( 含 灌浆 层 ) 的 
变形 ,以 及 锚 杆 与 岩 体 间 的 相互 作用 5 。 通 常 采 
用 拉 拔 实验 来 研究 锚 杆 与 岩 体 间 的 锚固 机 制 、 锚 杆 
轴 力 变化 和 界面 切 应 力 分 布 规律 等 。 

国内 外 学 者 围绕 锚 杆 拉 拔 的 理论 推导 和 实验 、 
数值 模拟 等 方面 开展 了 较 多 研究 。 如 李鹏 飞 等 '” 
基 扩 四 线性 黏 结 滑 移 界面 模型 ,对 拉 拔 力作 用 下 锚 
样 虹 面 的 渐进 失效 机 理 以 及 界面 的 切 应 力 分 布 规律 
进 伍 了 理论 分 析 ; 荣 冠 等 “通过 实验 研究 了 以 螺纹 
钢 秋 圆 钢 为 杆 体 的 全 长 秋 结 砂浆 锚 杆 的 锚固 机 理 和 
DERIER ADEE 对 锚 杆 和 砂浆 界面 的 黏 接 滑 移 
关 虹 进行 了 实验 研究 ,并 推导 了 界面 在 弹性 .塑性 软 
人 加 及 开裂 滑 移 阶段 的 力学 模型 ;Martin 4$ 4 Hh 
YET 锚 杆 拉 拔 端 位 移 和 销 杆 - 岩 体 界面 本 
MEG AUR A38 fi E 够 预测 拉 拔 实验 中 全 灌浆 钳 
FR 性 能 的 解析 方法 ; Ivanovic 457^ 基于 销 杆 与 
岩 硒 的 界面 作用 模型 ,对 错 杆 拉 氢 脱 茜 破坏 行为 进 
行家 研究 ,并 重点 探究 了 界面 届 服 和 软化 的 影响 。 
MEUS R FLAC” 模 拟 了 错 杆 拉 拔 中 不 同 围 压 、 
TR Q 值 以 及 岩 体 不 连续 面条 件 下 锚 杆 与 岩 体 间 的 
相 剖 作用 。 车 纳 等 ("采用 离散 元 软件 PFC 模拟 了 
绿 片 岩 内 的 锚 杆 拉 拔 过 程 ,对 不 同 锚固 长 度 、 围 压条 
件 下 锚 杆 的 破坏 机 理 进 行 了 研究 。 

在 数值 模拟 方法 中 , 非 连续 变形 分 析 ( discontin- 
uous deformation analysis , DDA ) 方 法 能 够 充分 考虑 岩 
体 的 不 连续 特性 ,在 岩 体 的 大 位 移 和 大 变形 分 析 方 
面具 有 独特 优势 ,已 成 为 岩石 力学 问题 分 析 的 重要 
数值 方法 31。 对 于 锚 杆 拉 拨 的 DDA 数值 模拟 ， 
Nie 等 4” 根据 拉 拨 过 程 中 灌浆 错 杆 与 岩 体 间 界 面 
的 力学 行为 ,推导 了 锚 杆 - 岩 体 耦合 作用 算法 并 集成 
到 DDA 当中 ,模拟 研究 了 不 同 端 部 条 件 和 界面 刚度 
对 锚固 效果 的 影响 ; Ma A P Hp EE HISP TR s 
滑 移 关系 简化 为 三 线性 模型 ,通过 DDA 模拟 锚 杆 拉 
拔 , 研 究 了 界面 切 应 力 的 变化 特征 和 界面 脱 务 ; 王 文 
xU" fg DDA 中 加 入 了 位 移 收敛 判断 准则 ,考虑 了 


锚 杆 的 塑性 变形 ,通过 锚 杆 拉 拔 模拟 验证 了 算法 的 
正确 性 ; Yokota 等 .采用 DDA 模拟 锚 杆 拉 拔 ,研究 
了 锚 杆 - 岩 体 界面 的 力学 行为 ,探究 了 锚 杆 肋 面 角 
度 、 肋 面 间隔 以 及 围 压 等 参数 对 锚固 效果 的 影响 。 

上 述 研 究 均 在 一 定 程度 上 考虑 了 锚 杆 拉 拔 中 的 
锚 杆 - 岩 体 界面 耦合 作用 ,但 对 锚 杆 拉 拔 的 力学 作用 
机 制 的 探讨 仍 有 待 深入 。 本 研究 在 DDA 中 引入 考 
虑 滑 移 脱 黏 的 锚 杆 - 岩 体 界面 耦合 作用 算法 ,通过 
DDA 模拟 系统 分 析 拉 拔 过 程 中 锚 杆 与 岩 体 间 的 界 
面 切 应 力 分 布 规律 与 界面 的 滑 移 脱 香 过 程 , 错 杆 轴 
力 的 变化 规律 ,以 及 他 们 与 错 杆 的 拉 拨 力 - 拉 拔 位 移 
曲线 间 的 关系 ,从 而 深入 揭示 锚 杆 拉 拔 的 力学 作用 
过 程 和 机 制 。 


1 锚 杆 拉 拔 的 DDA 模拟 方法 


图 1 为 锚 杆 拉 拔 DDA 模拟 模型 示意 图 ,在 此 基 
础 上 建立 锚 杆 - 岩 体 界面 耦合 作用 算法 。 锚 杆 被 n 
po: 点 划分 成 锚 杆 单元 ,第 1 个 和 第 n 个 锚 杆 
别 位 于 锚 杆 拉 拔 端 和 和 骨 入 端 ,并 以 拉 拔 端 为 
E x 轴 正 方向 建立 锚 杆 局 部 坐 
标 系 。 在 锚 杆 节点 相同 位 置 设置 岩石 节点 , 且 销 杆 
节点 和 岩石 节点 编号 相同 。 宕 体 划分 为 三 角形 
DDA 子 块 体 单元 ,以 模拟 岩 体 的 变形 。 由 于 在 
端 部 的 轴 向 拉 拔 作用 下 , 锚 杆 的 弯曲 变形 可 以 忽略 ， 
因此 对 于 锚 杆 ， aro di 
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NY AA 


错 村 NR 

, 图 1 锁 杆 拉 拔 DDA pp 

Fig.1 Illustration of DDA model to simulate rock bolt pullout 

在 锚 杆 拉 拔 实验 中 , 随 着 拉 拔 力 的 增加 , 锚 杆 与 

岩 体 间 界 面 的 切 向 变形 逐渐 增 大 , 发 生 剪 切 滑 移 并 

致 脱 黏 。 为 模拟 此 现象 ,引入 锚 杆 与 岩 体 界面 切 向 

耦合 作用 的 多 线性 务 结 滑 移 界面 模型 …” ,如 图 2 
所 示 。 


e——9 Hity 


ibi pett | 
到 2 多 线性 藉 结 滑 移 界面 模型 
Fig.2 Multilinear bond-slip interface models 

图 2 中 6 为 锚 杆 节点 与 岩 体 节点 的 切 向 相对 位 
移 ,7 为 节点 相 邻 锚 杆 单元 与 岩 体 间 的 界面 切 应 力 ， 
切 应 力 与 节点 相对 位 移 呈 多 阶段 线性 关系 。 在 三 线 
性 界面 模型 中 ,04 AB 和 BC 分 别 为 弹性 阶段 ,塑性 
软化 阶段 和 脱 莫 阶段。 在 四 线性 界面 模型 中 ,04、 
AB BC 和 CD 分 别 为 弹性 阶段 、 塑 性 剪 胀 阶段 、 塑 性 
软化 阶段 和 脱 莫 阶段 。 对 于 第 i 线性 段 , 定 义 模型 
参数 k; 和 c; 分 别 为 
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@@ 错 杆 拉 拔 过 程 中 , 锚 杆 单元 (长 度 为 dx) 的 受 力 
如 图 3 所 示 。 锚 杆 单元 与 岩 体 间 的 切 向 力 F., 和 错 
标的 轴 向 应 力 ou , 间 的 关系 为 

N (a. +do -or ) .A,= —F,,- dx (2a) 
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式 申 A, 为 锚 杆 横 截 面积 。 
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(b) 锚 杆 单元 受 力 


(a) "ii EE C HT 


图 3 锚 杆 模 截面 及 锚 杆 单元 受 力 示 意 医 


Fig.3 Illustrations of rock bolt cross-section 


and force analysis of bolt element 
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Ri pem 


AE, . [i A,E, AE, 
Rl i bo T i-1 PT i+l 
X 三 过 X =w N =N 
i-1 isl AE 
X 一 区 i pb isl 
; —sD:k()|-U - eU 
i- it i-1 it 
X -x X -x 


-nD-k(i)- Ui + 


“mTD:c(i) (3) 


ChinaXiv 合 作 期 刊 


sp LE, 为 锚 杆 弹性 模 量 ;x 为 节点 i 的 坐标 ; UnU, 
分 别 为 销 杆 节点 i 和 岩石 节点 i 沿 销 杆 轴 向 的 位 移 ; 
D 为 锚 杆 直径 。 对 于 直接 施加 拉 拔 力 下 的 第 1 个 节 
点 ,上 式 改 写 为 


A E po A E I 
uot mao) Meg 
X —X X —X 

x! 


2 
2 enD- c(1) € F 


(4) 
在 DDA HMF, MJI F Bi DDA 时 步 线性 增 
加 。 兰 石 的 节点 位 移 U, 采用 上 一 时 步 的 DDA 计算 
结果 ,进而 由 式 (3) 和 式 (4) 建立 锚 杆 节点 系统 平衡 
方程 组 ,并 采用 高 斯 消 元 法 求解 更 新 锚 杆 节点 位 移 
邦 。 当 锚 杆 节点 与 岩石 节点 的 相对 位 移 进入 脱 黏 阶 
Bt ,为 避免 方程 组 的 奇异 ( 脱 黏 阶段 相同 的 界面 切 应 
力 / 节 点 力 对 应 的 节点 相对 位 移 不 唯一 ) ,该 节点 的 
计算 不 再 进入 节点 系统 平衡 方程 , 即 节点 平衡 方程 
个 数 减 少 。 对 于 脱 竹 的 节点 ,以 相 邻 未 脐 秋 节点 为 
基准 ,通过 锚 杆 单元 伸 长 量 的 计算 来 更 新 其 位 移 。 
锚 杆 单元 的 伸 长 量 计算 如 下 。 


2 1 
IX. aD: k(1) A 


IY (5) 
RPF, HAMI; HATKE, IDET 
节点 与 岩石 节点 间 的 相对 位 移 , h E REREN 
界面 模型 插值 得 到 单元 界面 切 应 力 ,并 等 效 为 节点 
力 计 入 节点 系统 平衡 方程 组 的 玉 中 。 


2 ST REB) DDA 模拟 分 析 


首先 针对 文献 中 的 2 组 拉 拔 实 验 ,在 文献 中 的 
多 线性 界面 模型 参数 的 基础 上 ,基于 DDA 模拟 对 界 
面 模 型 参数 进行 修正 ,并 对 锚 杆 拉 拔 力 - 拉 拨 位 移 、 
轴 力 及 界面 切 应 力 进 行 模拟 和 分 析 。 此 外 还 针对 锚 
杆 的 弹性 模 量 进行 了 参数 分 析 。 


2.1 界面 模型 参数 修正 及 拉 拔 力 - 拉 拔 位 移 分 析 


刘波 等 ”进行 了 锚 杆 拉 拔 的 现场 实验 ,其 中 锚 杆 
长 度 工 =3 m, H1% D =42 mm, 弹性 模 量 E, =210 GPa; 
围 岩 材 料 为 细 石 混凝土 ,密度 p =2.1 x10 kg/m ,弹性 
Eit E, 225.6 GPa, 泊 松 比 凡 =0.26。 文 献 [ 13 ] 还 
在 实验 基础 上 开展 理论 推导 分 析 , 得 到 锚 杆 - 岩 体 界 
面 的 三 线性 界面 模型 参数 。 针 对 相同 实验 ,李鹏 飞 
等 则 经 过 理论 推导 分 析 得 到 了 锚 杆 - 宕 体 界面 的 
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在 明显 差异 。 因 此 ,本 研究 在 以 上 三 线性 和 四 线性 


四 线性 界面 模型 参数 , 均 列 于 表 1 。 

图 4 给 出 了 文献 [13] 及 文献 [11] 中 销 杆 拉 拔 
力 - 拉 拔 位 移 实 验 曲线 ,以 及 与 表 1 中 文献 [13 ] 及 文 
献 [11] 三 线性 及 四 线性 模型 参数 所 对 应 的 理论 模 
型 拉 拔 力 - 拉 拔 位 移 曲 线 。 可 以 看 出 ,理论 和 实验 所 
得 拉 拔 力 - 拉 拔 位 移 曲 线 存 在 较 大 差异 。 当 采用 文 
献 [13] 和 文献 [11] 中 的 模型 参数 时 ,通过 DDA 模 
拟 得 到 的 拉 拔 力 - 拉 拔 位 移 曲 线 与 实验 结果 同样 存 


界面 模型 参数 的 基础 上 ,在 DDA 模拟 中 对 模型 参数 
进行 不 断 调整 ,通过 DDA 模拟 得 到 的 拉 拔 力 - 拉 氢 
位 移 曲线 与 实验 曲线 进行 比较 ,分 别 得 到 修正 后 的 
三 线性 和 四 线性 界面 模型 参数 ( 见 表 1)。 从 图 4 可 
以 看 出 ,在 DDA 模拟 中 ,采用 修正 后 的 三 线性 和 四 
线性 界面 模型 均 能 更 好 地 吻合 实验 结果 ,并 且 四 线 
性 界面 模型 得 到 的 结果 与 实验 结果 贴 合 更 好 。 


表 1 多 线性 黏 结 滑 移 界面 模型 参数 


Tab.1 Parameters of the multi-linear bond-slip interface models 


模型 类 型 8,/mm v, / MPa 8,/mm 7,/MPa 63/mm 73/MPa 
文献 [13] 的 三 线性 模型 0.21 3.84 0.38 1.15 

修正 后 的 三 线性 模型 0.21 3.84 0.38 0.90 

CR C11] 的 四 线性 模型 0.10 2.11 0.38 2.45 0. 82 0.17 
d 

irs 四 线性 模型 0.10 2.11 0.25 2.45 0.60 0.90 
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(a) 三 线性 界面 模型 
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(b) 四 线性 界面 模型 
图 4 针对 文献 [13] 实 验 的 拉 拔 力 - 拉 拔 位 移 曲 线 对 比 


Fig.4 Pulling force-displacement curve comparisons regarding 


the referred experiment ^ 


在 图 4(a) 中 ,可 以 看 出 ,采用 修正 后 的 三 线性 
界面 模型 参数 得 到 的 拉 拔 力 - 拉 拔 位 移 曲线 可 较 明 
显 的 分 为 3 个 特征 段 。 第 I 阶段 (04 段 ) 为 弹性 阶 


段 , 曲 线 呈 线性 上 升 趋势 ,这 时 DDA 模拟 中 锚 杆 上 
所 有 节点 均 处 于 三 线性 界面 模型 的 弹性 阶段 ;第 工 
阶段 (4B8 段 ) 为 弹性 -塑性 软化 混合 阶段 ,这 时 靠近 
拉 拔 端的 节点 逐步 进入 三 线性 界面 模型 的 塑性 软化 
阶段 ;第 五 阶段 (BC 段 ) 为 弹性 -塑性 软化 - 滑 移 混 合 
阶段 ,这 时 靠近 拉 拔 端的 节点 逐步 进入 脐 黏 阶段 。 
在 图 4(p) 中 ,4.8 和 C 分别 代表 采用 修正 后 的 四 线 
性 界面 模型 时 ,靠近 锚 杆 拉 拔 端 节点 开始 进入 界面 
模型 的 塑性 剪 胀 阶段 .塑性 软化 阶段 和 脱 黏 阶 段 的 
时 刻 。 可 以 看 出 ,采用 四 线性 界面 模型 得 到 的 曲线 
更 加 平 请 ,其 阶段 性 不 明显 。 


2.2 锚 杆 界面 切 应 力 分 析 


针对 文献 [ 13 ] 的 拉 拔 实验 ,在 不 同 的 拉 拔 力 条 
件 下 ,分 别 采用 修正 三 、 四 线性 界面 模型 ,通过 DDA 
模拟 得 到 的 锚 杆 界面 切 应 力 曲 线 如 图 5 所 示 。 可 以 
看 出 , 锚 杆 界面 切 应 力 沿 锚 杆 轴线 的 分 布 形态 与 所 
采用 的 多 线性 界面 模型 形态 相似 , 且 随 着 拉 拔 力 的 
增加 , 销 杆 界面 切 应 力 峰 值 由 拉 氢 端 逐渐 向 区 入 端 
转移 。 就 销 杆 上 同一 位 置 而 言 , 随 着 拉 拔 力 的 施加 ， 
锚 杆 单元 /节点 移 处 于 弹性 阶段 , 当 采 用 四 线性 界面 
模型 时 ,随后 则 会 进入 塑性 剪 胀 阶段 , 直至 切 应 力 达 
到 峰值 ; 当 拉 拔 力 进 一 步 增 大 , 锚 杆 单元 /节点 随 之 
进入 塑性 软化 阶段 ,界面 切 应 力 开始 下 降 , 随 后 由 于 
界面 脱 香 的 发 生 ,界面 切 应 力 降低 至 恒定 残余 应 力 
而 不 再 变化 。 
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2/35 错 杆 轴 力 分 析 


< 一 在 不 同 的 拉 拔 力 条 件 下 ,分 别 采 用 修正 三 、 四 线 
VE iE ES ag sb DDA 模拟 得 到 的 锚 杆 轴 力 曲线 与 
IFRI 13 ] 的 实验 及 理论 模型 结果 对 比如 图 6 所 示 。 
可 以 看 出 , 锚 杆 的 轴 力 随 着 嵌入 深度 增加 而 逐渐 减 
小 5 址 轴 力 曲线 形态 特征 会 随 着 拉 拔 力 的 变化 而 改 
变 。 图 中 的 软化 分 界 点 为 锚 杆 上 最 大 界面 切 应 力 所 
在 位 置 ,开裂 分 界 点 为 距离 拉 拔 端 距 离 最 远 的 处 于 
脱 黏 阶段 的 节点 位 置 。 因 此 ,开裂 分 界 点 往 左 与 锚 
杆 拉 拨 端 之 间 的 错 杆 单元 处 于 脱 莫 阶段 ,开裂 分 界 
点 往 右 至 软化 分 界 点 之 间 的 锚 杆 单元 处 于 塑性 软化 
阶段 ,而 软化 分 界 点 往 右 至 圣人 端的 锚 杆 单元 处 于 
弹性 或 塑性 剪 胀 阶段 。 

从 图 6 可 以 看 出 ,在 不 同 的 拉 拔 力 条 件 下 ,DDA 
模拟 结果 与 实验 结果 均 吻 合 较 好 。 当 拉 拔 力 下 = 
160 kN 时 , 销 杆 上 所 有 节点 均 处 于 弹性 阶段 ,节点 
力 随 嵌入 深度 持续 减 小 , 轴 力 从 拉 拔 端 到 矢 人 端 逐 
渐 下 降 , 曲 线 总 体 呈 现下 四 形态 。 当 下 =250 kN 时 ， 
轴 力 曲线 在 软化 分 界 点 左 侧 旦 上 凸 形态 ,这 是 因为 
该 范围 内 节点 处 于 塑性 软化 阶段 ,上 且 越 靠近 拉 拔 端 ， 
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锁 杆 节点 力 越 小 。 曲 线 在 软化 分 界 点 右 侧 时 下 四 形 
态 ,这 是 因为 该 部 分 节点 处 于 弹性 或 塑性 剪 胀 阶段 ， 
越 靠近 软化 分 界 点 , 锚 杆 节点 力 越 大 。 与 7 = 
250 kN 类 似 , 当 =300 kN 时 , 轴 力 曲线 在 软化 分 界 
点 左 侧 呈 凸 形态 , 右 侧 呈 下 四 形态 。 开 裂 分 界 点 
左 侧 的 节点 处 于 脱 黏 阶段 ,这 时 锚 杆 单元 界面 切 应 
力 为 恒定 的 残余 应 力 。 
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Fig.6 Axial force variation curves of the bolt 
2.4 锚 杆 弹性 模 量 对 模拟 结果 的 影响 


上 文 锚 杆 拉 拔 力 - 拉 拔 位 移 、 界 面 切 应 力 和 轴 力 
的 分 析 表 明 ,采用 本 研究 所 发 展 的 销 杆 拉 拔 DDA 模 
拟 方法 ,基于 修正 的 多 线性 界面 模型 ,能够 较 好 地 再 
现 拉 拔 实验 力学 作用 过 程 ,并 在 此 基础 上 揭示 锚 杆 
拉 拔 的 力学 机 制 。 以 下 在 文献 [13 ] 的 实验 基础 上 ， 
采用 修正 后 的 四 线性 界面 模型 ,对 不 同 锚 杆 弹性 模 
量 条 件 下 的 锚 杆 拉 拔 进行 模拟 ,分 析 锚 杆 弹性 模 量 
对 拉 拔 力学 响应 的 影响 。 图 7 为 DDA 模拟 得 到 的 
锚 杆 拉 拔 力 - 拉 拔 位 移 、 界 面 切 应 力 和 轴 力 结果 ,其 
中 销 杆 弹性 模 量 分 别 取 值 为 0. 9E, , 0. 95E, 
1. 00E,,1.05E, fI 1. 1E,CE, 2210 GPa), 
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图 7 不 同 锚 杆 弹性 模 量 条 件 下 的 模拟 结 
Fig.7 Simulation results under different bolt modulus conditions 

从 图 7(a) 可 以 看 出 , 随 着 锚 杆 弹性 模 量 的 增 
加 ,相同 拉 拔 力 条 件 下 的 拉 拔 位 移 减 小 , 即 锚 杆 的 轴 
向 伸 长 量 也 减 小 ; 同时 , 随 着 拉 拔 力 的 增 大 ,由 于 靠 
近 拉 拔 端 的 节点 逐渐 进入 塑性 软化 和 脱 黏 阶 段 , 减 
小 了 岩 体 对 错 杆 的 切 向 约束 作用 ,弹性 模 量 对 拉 拔 
位 移 的 影响 也 越 来 越 大 。 此 外 ,在 相同 的 拉 拔 力 
(F =300 kN) 条 件 下 ,由 于 较 大 的 弹性 模 量 导 致 更 
小 的 锚 杆 伸 长 量 , 减 小 了 锚 杆 与 岩 体 间 的 相对 位 移 ， 
从 而 使 得 锚 杆 上 靠近 拉 拔 端 达 到 塑性 软化 和 脱 黏 阶 
段 的 节点 数量 相对 较 少 ,对 应 的 节点 力 更 大 。 因 此 ， 
如 图 7(b) 和 图 7(c) 所 示 , 随 着 锚 杆 弹性 模 量 的 增 
加 ,界面 切 应 力 峰 值 更 加 靠近 拉 拔 端 , 且 沿 着 锚 杆 拉 
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拔 端 向 垦 入 端的 方向 , 销 杆 轴 力 下 降 的 速率 更 大 。 


3 结 it 


本 研究 在 DDA 中 引入 了 锚 杆 - 岩 体 界面 而 合作 
用 算法 ,采用 多 线性 黏 结 滑 移 界 面 模 型 模拟 锚 杆 的 
消 移 脱 黏 ,针对 实际 拉 拔 实验 开展 对 照 模拟 ,研究 了 
锚 杆 拉 拔 的 力学 作用 过 程 和 机 制 。 所 得 结论 如 下 。 

1) 基 于 修正 参数 的 多 线性 界面 模型 ,模拟 得 到 
的 销 杆 拉 拔 力 - 拉 拔 位 移 、 界 面 切 应 力 和 轴 力 的 变化 
和 分 布 规律 与 实验 结果 吻合 良好 ,验证 了 所 提 算 法 
的 正确 性 与 适用 性 。 

2) 采 用 三 线性 界面 模型 时 , 拉 拔 力 - 拉 拔 位 移 曲 
线 呈 现 较为 明显 的 3 阶段 特征 ,而 采用 四 线性 界面 
模型 时 , 拉 拔 力 - 拉 拔 位 移 曲线 更 加 平滑 , 且 与 实验 
结果 吻合 程度 更 高 。 

3) 锚 杆 界面 切 应 力 的 分 布 形态 与 多 线性 界面 模 
型 形态 相似 , 沿 租 入 端 向 拉 拔 端 逐渐 经 历 弹性 和 塑 
性 剪 胀 、 塑 性 软化 及 脱 黏 阶段 。 随 着 拉 拔 力 的 增加 ， 
界面 切 应 力 峰 值 逐渐 由 拉 拔 端 向 能 和 人 端 转 移 。 

4) 在 不 同 拉 拔 力 条 件 下 , 锚 杆 轴 力 曲线 上 的 软 
化 分 界 点 和 开裂 分 界 点 位 置 与 实验 结果 基本 一 致 。 
拉 拔 力 较 小 时 , 锚 杆 界面 处 于 弹性 或 塑性 剪 胀 阶段 ， 
轴 力 曲线 呈 下 四 形态 ; 随 着 拉 拨 力 的 增 大 ,靠近 锚 杆 
拉 拔 端的 节点 逐渐 进入 塑性 软化 或 脱 黏 阶段 ,对 应 
部 分 的 轴 力 曲线 呈 上 凸 形态 。 

5 ) 随 着 销 杆 弹性 模 量 的 增加 ,由 于 锚 杆 轴 疝 伸 
长 变形 量 的 减 小 ,相同 的 拉 拔 力 将 导致 更 小 的 拉 拔 
位 移 , 且 界 面 切 应 力 峰 值 更 靠近 拉 拔 端 , 锚 杆 轴 力 治 
拉 拔 端 向 钥 入 端 下 降 更 快 。 

本 研究 锚 杆 - 岩 体 界面 耦合 作用 算法 的 实现 ,也 
为 岩 体 锚固 的 DDA 模拟 分 析 创 造 了 条 件 。 
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